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TWO BONDED HALF PLANES WITH A 
CRACK GOING THROUGH THE INTERFACE* 


BY 

F. Erdogan and V. Biricikoglu 
Lehigh University/ Bethlehem, Pennsylvania 


ABSTRACT 

The plane problem of two bonded elastic half planes con- 
taining a finite crack perpendicular to and going through the 
interface is considered. The problem is formulated as a system 
of singular integral equations with generalized Cauchy kernels. 
Even though the system has three irregular points, it is shown 
that the unknown functions are algebraically related at the 
irregular point on the interface and the integral equations 
can be solved by a method developed previously. The system of 
integral equations is shown to yield the same characteristic 
equation as that for two bonded quarter planes in the general 
case of the through crack, and the characteristic equation for 
a crack tip terminating at the interface in the special case. 
The numerical results given in the paper include the stress 
intensity factors at the crack tips, the normal and shear com- 
ponents of the stress intensity factors at the singular point 
on the interface, and the crack surface displacements. 
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1. INTRODUCTION 


In considering the fracture of composite materials, one 
may approach - the probl em fronT'two ~ diff erelft'p oi n't s”of view. 

In the first approach the primary interest is in studying and 
in estimating the "bulk strength" of the given structure under 
a known system of external loads and environmental conditions. 
In this type of studies it is usually assumed that the existing 
imperfections such as voids, cracks, and inclusions are ran- 
domly distributed throughout the composite medium and the mate- 
rial is statistically homogeneous. Thus, the very nature of 
the problem requires that some kind of a statistical strength 
theory be used as a guide in the investigations. 

In the second approach to studying the fracture of com- 
posites, one is basically interested in the initiation of frac- 
ture from the "localized" imperfections which are known (or 
assumed) to exist in the material. In this type of studies it 
is usually assumed that the composite medium consists of per- 
fectly bonded elastic components and the localized imperfection 
may be idealized as a plane crack or as a flat elastic inclu- 
sion. The problem of a finite crack lying parallel to or at 
a bimaterial interface in composites with various geometries 
was discussed in [1-5]. The similar problem for a flat inclu- 
sion was studied in [6]. An up-to-day review of the general 
fracture problems in composite materials and a summary of some 
of the known results may be found in [7]. The distinguishing 
feature of the solutions given in [1-7] as well as the other 
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known solutions which appeared in literature within the past 
decade is that the strength of the stress singularity at the 
imperfection front is -1/2 and, in the case of a crack, the 
quasi-static stress state in the neighborhood of the crack 
front remains autonomous as the crack propagates. That is, 
aside from a slight change in a multiplicative constant known 
as the stress intensity factor, the asymptotic nature of the 
stress state in the vicinity of the crack front remains 
unchanged . 

On the other hand, when the crack front terminates at a 
bi-material interface in the composite medium the strength of 
the stress singularity is no lorvger -1/2 and the angular dis- 
tribution of stresses differs considerably from that of a crack 
tip imbedded in a homogeneous medium. After reaching the inter- 
face, further propagation of the crack may be in the form of 
(a) a cleavage crack into the second medium (Figure 1), (b) a 
debonding crack along the interface, or (c) a "reflected" crack 
back into the first medium. The basic problem of a finite 
crack terminating at a bi-material interface was discussed in 
some detail in a previous paper [8]. [8] also presents a ten- 
tative fracture criterion which may be used in predicting the 
mode of the fracture propagation (such as (a)-(c), mentioned 
above) as well as the level of the external loads initiating 
the fracture provided the fracture strengths of the two adjoin- 
ing materials and of the interface are known. 

In this paper we consider the problem of a crack going 
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through the interface into the second medium. From the analyt- 
- ica1_v_ie.w.po,int— the unusual feature of the p rob lem is that it 
has three irregular points (two crack tips and the intersection 
of the crack with the interface) at which the stress state is 
generally singular. One of the objectives of the paper is the 
development of the necessary theoretical and numerical tech- 
niques for handling problems of this nature. 

2. THE INTEGRAL EQUATIONS 

Consider the plane elasticity problem shown in Figure 1. 

In this paper only the symmetric problem will be considered; 
that is, it is assumed that in the composite medium without 
the crack and subjected to the same external loads as the 
cracked medium, the planes (r > 0, 0 = ir) and (r > 0, 0 = 0) 
are free from shear and are acted upon by normal stresses p-|(r) 
and P 2 ( r ) » respectively. Thus, through a proper superposition, 
the singular part of the problem may be reduced to that shown 
in Figure 1, where the sel f -equ i 1 i brati ng crack surface trac- 
tions - p 1 ( r ) and - p 2 ( ^ ) are the only external loads. In the 
special case of uniform loading p^ and p^ are constant and the 
continuity conditions require that ^ 

p l p 2 

for plane stress, 

t l t 2 

1 - v, 2 1 - \>j- 

— f — - — p, = — f — — p, , for plane strain. (l.a,b) 

t l 1 k 2 £ 

To solve the problem shown in Figure 1 we will first 
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consider the problem described in Figure 2 where each half 
_plane contains- a finite crack perpendicular to the interface, 
and then let the distances a^ and a 2 go to zero. 

Because of symmetry, it is sufficient to consider the 


problem in the half 
related elasticity 
boundary conditions 

plane 0 < r < <», 
problem has to be 

0 < 0 < tt (Figure 
solved under the 

2). The 
f ol 1 owi ng 

T 1 re ( r ,1T ) = 0 

8 

V 
JL 

VI 
O 



T le e(r,Tr-°) = 

- P] (r) , a 1 < 

r < b 1 , 


U 1 6 ( r > 7r ) = 0 » 

0 < r < a-j , 

b.| < r < °°; 

(2 . a-c ) 

U 16 = u 20 * 

u lr = u 2r * T 100 

t 200 * 


T lre = T 2r6 * 

CD 

II 

=* 

ro 

o 

< r < ; 

( 3 . a-d ) 

T 2re (r,°) = 0 

8 

V 

VI 

o 



T 2ee^ r,+0 ^ 

P 2 (r) » a 2 < 

r < b 2 . 


u 2 0 ( r » 0 ) * 0 » 

0 < r < a 2 » 

b 2 < r < . 

(4. a-c) 

Let the unknown functions f-| and 

f 2 be defined as 

f ol 1 ows : 

V r > - - 2 If 

u l0 (r,ir-O) , 



f 2 (r) ’ 2 If u 

26 (r,+°) . 


(5.a,b) 


From (2.c), (4.c) and (5) it follows that 



f-j(r) = 0 , 0 < 

f 2 (r) = 0 » 0 < 

b l 

/ f -i ( r)dr = 0 , 
a l 


r < 


u 2 

/ 

a 2 


a-j , b-j < r < 
a 2 , b 2 < r < «>, 

f 2 (r)dr = 0 . 


( 6 .a-d) 


Using Mellin transforms and following the procedure out- 
lined in [8], after somewhat routine manipulations the problem 
can be reduced to the following system of integral equations 
for the functions f-j and f 2 : 


i b l Ms) , b l 

? I -W ds + ? I k 11 (s)ds 


l 2 1+K i 

+ n f k 12 (r,s)f 2 (s)ds = - — P 1 ( r ) 

a o 1 


a, < r < b, , 


, b 2 f ? (s) , b l 

7 / “sTF- ds + ? / ^(r.sjf^sjds 


K 2 


l 2 1+K ? 

+ ^ / k 22 (r,s)f 2 (s)ds = 

a 2 


where the kernels k . • are given by 

' J 


P 0 (r) , a 2 < r < b 2 , 
( 7 . a ,b) 


k 1 1 ( r , s ) = l 


3 c lk r 


k-1 


k=l (s+r) 


F 


k 2 2 ( r » s ) = J 


3 c 2k r 


k-1 


k=l (s+r) 


F 


k 12 (r»s) = 


? d lk r 


k-1 


(s+r) 


2 d 2k rk_1 

k 2 i ( r » s ) = | ( s+r )k ; (8 * a " d) 
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1 mjfl+Kj) 3 ( 1 -m 1 ) 

C 1 1 2 2((n-| +k 2 ) 2 ( 1 +m-| ) * 


12 


6(1-1,) 

1 +m, ic ^ 


13 


4(1 -m, ) 
T+nTpioj 


•, m 2 (l+ic 2 ) 3(l-m 2 ) 

'21 F 2(m 2 +Kj ) ” 2( 1 +m 2 ic 2 ) * 


'22 


6(1 -m 2 ) 

1 +m 2 < 2 


23 


4(1 -m 2 ) 
1 +m 2 K 2 


3( 1 +K-J ) 1 +K-| 

*11 “ 2(m 2 +K 1 ) ” 2(l+m 2 K 2 ) * 


'12 


1 +K-j 

1 +m 2 K 2 


1 +K, 


m 2 +K, ’ 


3(1+k 2 ) 


1 + K, 


21 2 (m.| +< 2 ) 2 ( 1 +m^ k, ) 


1 +K, 


1 + K, 


22 ~ 1 +ni-j K-j m^ +< 2 


m l = y 2^ y l * 


m 2 p 1 /y 2 * 


(9) 


In (9) y ^ and y 2 are the shear moduli, = 3-4v^ for plane 
strain and k.. = (3-v, )/(l+v. ) for the generalized plane stress, 
(i = 1,2), v.. being the Poisson's ratio. 


From (8) it is seen that for a^ > 0, a 2 > 0 the kernels 
k..(r,s), (i,j = 1,2) are bounded and continuous. Hence (7) 

I J 

is a simple system of singular integral equations with an index 
k = 1 and the additional conditions (6.c and d) which are re- 
quired for a unique solution. These equations can be solved 
for the unknown functions f-j and f 2 in a standard manner. A 
simple and highly effective numerical method has been described 
in [9]. After obtaining f^ and f 2 , all the desired field 

1 

quantities may be expressed as and evaluated from the definite 
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integrals with kernels corresponding to the particular field 
quantity and the density functions f , , (i = 1 , 2 ) (see, for 
example, [ 8 ] for the case of one crack). From the practical 
viewpoint of particular interest are the stress intensity fac- 
tors defined by* 


k ( a i ) 


r 2y, 

lim /2( a, -r) T 1 QB (r,Tr) = lim TT — ^(r-ajf, (r) , 
r+a, 1 100 r-a, ,+lc l 11 


Mb,) 


= lim / 2 (r-b, ) 

r+b, 


X] e0 (r .tt) 


2y, _ 

-lim TT 7 - / 2 (b, -r) f,(r) 
r+b, 1 K 1 ' 1 


_ 2 y ? 

k ( a g ) = lim / 2 (a 2 -r) T 2 e 0 (r,O) = lim / 2 ( r-a 2 ) f 2 (r), 

r~ > 32 2 

k(b„) = lim / 2 (r-b ? ) x ?flfl (r, 0 ) = -lim yr|- / 2 (b ? -r) f ? (r) 
c r+b 2 L ^ r^b 2 1 K 2 c * 

(lO.a-d) 

For the bonded planes Material 1 Aluminum (E, = 10^ psi , 
v, = 0.3), and Material 2 Epoxy (E 2 = 4.45*10^ psi, v 2 = 0.35), 
Tables 1 and 2 show some of the calculated plane strain results 
In these tables p, and p 2 refer to constant pressures applied 
to the crack surfaces and are related by (1). As a basis of 
comparison the special case of a homogeneous plane containing 
two collinear cracks may be mentioned. For y, = y 2 , v, = v 2> 
a, = a 2 = a, b, = b 2 = b, the stress intensity factors are 
given by 


* 


See Section 4 below. 
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Table 1. 


Stress Intensity Factors for a-j >0, a 2 > 0. 

Ma teri al T r ATumi numT— Materi al— 2 E-po-xy ; 

a 2 = 1 » ^2 = (^2“ a 2^^ ~ ^ ^ » 

c-j = ( b -j +a -j ) /2 = 1.5, A-j = (bj-a-|)/2 variable. 


£ 1 

k(a 1 ) 

k(b 1 ) 

k(a 2 ) 

k(b 2 ) 

P 1 * / *7 

Pl^T 

P 2 ^2 


0.05 

1 .0015 

1 .0014 

0.97417 

0.98165 

0.10 

1 .0023 

1 .0021 

0.97515 

0.98241 

0.15 

1 .0037 

1.0033 

0.97679 

0.98369 

0.25 

1 .0082 

1 .0069 

0.98212 

0.98775 

0.50 

1 .0326 

1 .0231 

1 .0086 

1 .0079 

0.75 

1 .0836 

1 .0492 

1 .0590 

1.0454 

1 .00 

1.1832 

1 .0869 

1.1474 

1.1089 

1 .25 

1 .4094 

1 .1437 

1.3118 

1.2210 

1 .40 

1 .8046 

1 .2012 

1 .5047 

1.3451 


Table 2. Stress Intensity Factors for a-| >0, a 2 > 0. 
Material 1: Aluminum, Material 2: Epoxy; 
a 2 = J^i = = 0*5. c, = ( b -j +a ^ ) / 2 variable. 


C 1 

k ( a 1 ) 

k(b 1 ) 

k ( a 2 ) 

k ( b 2 ) 

p r^T 

p i 

p 2 2 

P 2 ^ ^ 2 

1 .50 

1 .0326 

1 .0231 

1 .0086 

1 .0079 

1 .25 

1 .0494 

1 .0326 

1 .0187 

1.0143 

1.00 

1 .0844 

1.0499 

1.0344 

1 .0236 

0.75 

1.1849 

1 .0881 

1.0624 

1.0393 

0.60 

1.4114 

1.1454 

1 .0965 

1.0576 
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k(a) _ [b 2 E(m)/K(m) - a 2 ]/2 

p Q /T ( b-a ) /a ( b+a ) 

MM. = /? b 2 [1 - E(m)/K(m)] () , , b) 

p 0 /l (b-ajv'btb+a) 

U = ( b - a ) / 2 , m = 1 - a 2 /b 2 ) , 

where p Q is the constant pressure applied to the crack surfaces 
and K and E are the complete elliptic integrals of first and 
second kind, respectively. For example, for b = 2a which corre- 
sponds to the row = 0.5 in Table 1 (or to c^ =1.5 in Table 
2), we have k ( a ) = 1.01762 P Q /£> k ( b ) = 1.01249 P 0 /1F. 

As ->1.5 in Table 1 and- c-| -> 0.5 in Table 2, k(a-j) •> °° 
while the remaining three stress intensity factors approach 
some finite values. The reason for this is that in this limit- 
ing case a i = 0 and the singular behavior of the stresses around 
the crack tip is of the form R a with 0.5 < a < 1 whereas the 
definition of k(a-|) is based on R~ 0 ' 5 type singularity (see (10)) 
where R is a small distance from the crack tip a-j (see [8]). 

The results for another special case are shown in Figure 3 where 
~ b 2' This is the problem considered in [8] where only one 
of the adjoining materials contains a crack. Because of the 
change in the power a of the stress singularity when the crack 
terminates at the interface, again as a + 0 k(a) -> °° for 
u-| > and k ( a ) 0 Tor y i < y 2 » whereas k(b) in both cases 
remains bounded. 
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3. CRACK GOING THROUGH THE INTERFACE, a-| = 0 = a 2 

We now consider the main problem in which the crack goes 
through the interface, occupying both materials. In this case 
it is clear that the integral equations (7) are still valid, 
with a^ = 0 = & 2 ' Thus, the equations have only three irregular 
points ( b -j , tt ; 0; b^.O), meaning that one has to be more careful 
in examining the singular behavior of the solution. Unlike the 
problem considered in the previous section for which a-j > 0, 

&2 > 0» as seen from (8), here the kernels k^. are no longer 
bounded for all values of r and s in the closed intervals [ 0 , b ^ ] 
and CO , b 2 H . At the irregular point 0 they become unbounded, 
indicating that at r = 0 the behavior of f-(r), (i = 1,2), will 
be affected by k..(r,s). 

* J 


To examine the behavior of f^ and f£ around the irregular 
points, following [10] we assume that the unknown functions may 
be expressed as 


V s) 


iri 8 . 

9j(s) _ 9 j ( s ) e 

s a (b r s) 6j s“(s-b j ) 6j 


( 12 ) 


0 < Re (a) < 1 , 0 < Re( B • ) < 1 , (j = 1,2), 

J 


where a, B-j and B 2 are the powers of the singularity at the 
three irregular points, the function g- satisfies a Holder 

J 

condition in the closed interval 0 < s < b., (j = 1,2), and 

n 3i 

the function s (s-b.) J is any definite branch which varies 

J 

continuously on 0 < s < b i# (j = 1,2). Consider now the fol- 

J 

lowing sectionally holomorphic functions: 


13 


i Ms) , b j g, 

6 . ( z) = - f -J ds = - / — i 

s ._„ z 7 Hr— 


iri 3, 

(s)e J ds 


— 3~. 

s a (s-b.) J (s-z) 




.1,2),. 

(13) 


where the complex variable z = r + ip is defined in such a way 
that in each case the cut (p = 0, 0 < r < b . ) , (j = 1,2) lies 

vl 

along the positive real axis. Separating the leading terms at 
the end points of the cuts, (13) may be expressed as [10] 


g^oie 

g j ( z) = 


in a 


i W 

b'Jsm™ z “ b j as,n ” e 

vJ 


1 


j (z-bj) 0 


67 +G oj( z >- < 14 > 


where G . is bounded everywhere except possibly at the ends, 

J 

e j ^ , near which 


i G oj (z > I < 


h . l 


z - e jk 


jk 


, ( J ,k = 1,2) 


(15) 


e jl = °» P jl < Re (“) ; e j 2 b j* P j 2 < Re ^j)» 

and Pj k being real constants. 


Using now the Plemelj formula [10] 


1 J f i ( s ) , 

i/ ds “ f [G j (r) + G j" (r)] > ( 0 < r < bj ) (16) 

the terms with the Cauchy singularity in (7) may be expressed 
as 


14 



1 J f i (s) 9.(0) ! 9i(M ! 

? « ds ' cot ’° 7 ’ IT" cot ^ ^71 

J J 


+ G ij (r) • 


(0 < r < b jS j = 1 ,2), 


07) 


where the behavior of G..(r) in the closed interval [ 0 , b - ] is 

* J J 

similar to that of G q j (j = 1,2). The remaining terms in the 
integral equations (7) may be obtained directly from (14) 
through differentiation. From (7), (8) and (13) it is seen 
that in these terms z = -r which is outside the cut 0 < r < b., 

J 

j = 1,2. As r varies in the interval (0,b.), z will vary in 

J 

the interval ( - b - , 0 ) in which G.(z) is holomorphic, (j = 1,2). 

J J 

Thus, from (13) and (14) we obtain 


i b J Ms)ds g -: ( 0) 

* V ■ V- r > • -rj— r 

0 b. J simra r 

J 


1 a + ^ 2 j ( r ) ’ 


1 3 rf.(s)ds . ag.(O) , 

? 77^ = " r ^ G J ( " r) = “tJ 7 + G ^ (r) 4 

0 ' s n b i J simra r 

J 


i 7 r 2 f i (s)ds 2 ,2 a(a+l)g.(0) , 

• I } 3 ■ 5- br G j<- r > ■ B. J hi * G 4j( r > ■ 


(s+r) 


T~ 


dr 


2b. ^sinua r 


( 0 < r < b . , j = 1,2), ( 1 8 . a - c ) 

J 

where it may easily be shown that the behavior of G ^ j ( r ) in 
the closed interval [0 ,b (j = 1,2, k = 2,3,4), will again 

J 

be similar to that of G qj . given by (15). 


If we now substitute from (17), (18) and (8) into (7), 
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and multiply first both sides of each equation by r a and let 

8i 

r -*■ 0, then (7. a) by (b,-r) and let r -*■ b^ , and (7.b) by 
(b~ 2 - r ) ^ and let r b 2 , we obtain the following system of 
equations : 


( c o s ttcx + c 1 1 + ac 12 + ^- ” - +1 ) c i 3 ) — ^ 


9i (0) 


b^ ^sinira 


+ (d^ + ad 12 ) 


9,(0) 


= 0 , 


\>2 sinira 


( d 2 1 + ad 22 ^ & 


9,(0) 


b^ ^sirnra 


+ (cosira + c 2l + ac 22 + a ^ +1 ) c,,) — tt — - 


2 -23 1 e 2 

b 2 simra 


0 , 


9, (b, ) 

— cotTr3 n = 0 

. Ol I 


9 2 (b 2 ) 


b 


a 

2 


cotir3 2 


0 . 


(1 9.a-d) 


g.(0) and g.(b.), (j = 1>2), are non-zero constants which will 

J J J 

be shown to be related to the stress intensity factors. Thus, 
(19. c) and (19. d) give the following characteristic equations 
to determine 3-j and $ 2 : 

cotir$i = 0 , cotTrB 2 = 0 . (20. a, b) 

From (20) the values of 3, and 0 2 satisfying the condition 
0 < Re(3j) <1, (j = 1,2), are found to be 3-j = 3 2 = l/2,!which 
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is the well-known result. 


On the other hand, since a t 0 and a / 1, for (19. a) and 
(19. b) to give a non-zero solution for g.(0), the determinant 

J 

of the coefficients must vanish, giving the characteristic 
equation to determine a as follows: 

D(a) = (cosua + c^ + ac -| 2 + — ^ - ^ - + c i (costto. + C21 

aC 22 ^2 C 23 ^ ^ 1 2 ) ( d 2 1 ac *22 ^ = ® • 

( 21 ) 

From (19) it is clear that the constants g-|(0) and g 2 ( 0 ) are 
not independent and are related by 


9 2 (°) = " 9] (°) » /b 2^ b l 


cosira + c 


11 


+ ac 


12 


a( a+1 ) 

2 


13 


'll 


+ ad 


12 


( 22 ) 

As will be pointed out later, the use of (22) will be necessary 
in order to obtain a unique solution for the integral equa- 
tions (7). 


In the special case of homogeneous medium, m . 1 , 

J 

<1 -*■ k 2 , c. k -► 0 (j = 1,2; k = 1,2,3), d^ -> 1, d ]2 -► 0, 

d 2 i = 1. d 22 -*■ 0, the characteristic equation (20) becomes 
2 

cos -rra - 1 =0, giving a = 0 as the acceptable root and 
g 2 (0)//B^ = - g-|(0)//5y. Considering the definitions (5) and 
(12), this indicates that, as expected, at r = 0 the derivative 
of the crack surface displacement becomes finite and continuous. 
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In the problem of two bonded half planes containing a 
through crack perpendicular to the interface (Figure 1), the 
singular behavior of the solution at and around r = 0 must be 
the same as that for two bonded quarter planes (see, for exam- 
ple, [11] and [12]). Defining the constants 

K - ! = 2 ( m i - 1 ) , K 2 = m-| (1 +Kj ) - ( 1 "*"^2 ) * 

K 3 = m^(l+K-|) + (I+k^) » (23) 

after some manipulations ( 21 ) may be reduced to 

[(K r K 2 ) cos 2 p - K,(a-1 ) 2 ] 2 

+ Kj 2 sin 2 p cos 2 p- K , 2 ( a- 1 ) 2 - 0 . (24) 

(24) is identical to the characteristic equation found in [11] 
and [ 12 ] which indicates that the two solutions will have the 
same singular behavior*. 

In the special case of only one of the cracks terminating 
at the interface, i.e., for b-j > a^ = 0 , bg > ^2 > the 
kernels k^* ^21 anc * ^22 in will be bounded and will have 
no contribution to the singularities. For this problem the 
solution Is of the form 


As shown in [11] the weak power singularity 0 < Re(a) < 1 as 
assumed in this paper is not the only possible singularity for 
two quarter planes. Theoretically it is possible to have a 
logarithmic singularity as well as bounded stresses. However, 
if one studies the conditions in detail, it is not difficult 
to show that for all practical material combinations, the 
power singularity is the only possibility. 
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Vs) = 


9 


f 2 (s) = 


9! (s) 


a, 3 

s 1 (b 1 -s) 


1 


g 2 (s) 


Ot o 3 o 

(s-a 2 ) 2 ( b 2 - s) 2 


9 


(25. a, b) 

and, following the procedure outlined above, the characteristic 
equations for a. and 3^ (i = 1,2) are found to be 

C 0 STra-| + c -j -j + ac-| 2 + c 13 = 0 , 

cot7r3-| = 0 , cotuoig = 0 , cotir32 = 0 • (26.a-d) 


(26. a) is the equation found in [8] for the single crack and 
( 26 . b-d ) give 3-| = a 2 = 3 2 = 

4. STRESS INTENSITY FACTORS 


Using the results of the previous section it can be shown 
that the strength of the stress singularity at the crack tips, 
conventionally known as the stress intensity factor, is related 
to and can be evaluated from the density functions or the dis- 
placement derivatives f.(r), (j = 1,2), defined by (5). First, 

J 

consider the case of two cracks imbedded into the adjacent 

homogeneous materials (i.e., b. > a i > 0, j = 1,2) for which 

the functions f- and the related holomorphic functions G-(z) 

J ^ 

may be expressed as 


fj(s) 


9j (s) 


a • 3 . * 

(s-a . ) J( b ;-s) J 


J 


(27. a) 
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'W. 



1 


G j( z ) 


* l 


b j fj(s)ds 


S-2 


iria . 

9 ,i (a .i> e 


J 


3.: a • 

(bj-aj)— J simTa J .(z-aj)~^ 


hlhL 


1 


a a 3 .• 

( b j' a j) JsinTre j ( z ' b j) J 


+ 6 Q j(z) , ( j = 1 ,2) 


(27. b) 


where = 3j = 1/2. We now note that in equations (7) the 

left hand sides are the expressions for [ ( 1 +Kj )/2y . ]t j 0q for 

(0 < r < a., b. < r < °°) as well as for (a. < r < b) as indi- 
J J J J 

cated in (7). Thus, noting that the terms involving k. . are 

* J 

bounded, outside the cuts G j ( r ) directly gives Tjqq on the real 
axes, from which the stress intensity factors may be readily 
evaluated. For example, at the end r = b-j from (27) we obtain 

k(b,) = lim /2(r-b, ) T lflfl (r,ir) = lim /2(r-b, ) G,(r) 

1 r-*-b-| 1 Iay 1+K 1 r-*-b ] 1 1 


2p, 9j (b ] ) 

1+IC 1 /b ] -a 1 


2y 


1 


1 +K 


1 1m 
1 r-*-b 


1 


/2 ( b-j -r) f ] (r) . 

(28) 


The other expressions given in (10) follow from (27) in a simi- 
lar way which, in terms of g., become 

J 


k(b 2 ) 


2y 2 /2 g 2 ( b 2 ) 

1 +IC 2 /b 2 -a 2 


2y, /? g • ( a ■• ) 

k ( a i ) = J **— J- , (j = 1,2) (29. a, b) 

Consider now the second typical case, that is, a crack tip 
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terminating at the interface (e.g., b-j > a^ = 0, > ag > 0) . 

In this case the stress intensity factor will be defined as 


k i ( 0 ) = lim /l r 1 T 200 (r,O) , 


(30) 


where a-j is obtained from (26. a) and T 2 ee( r » 0 ) given by 
(7.b). Since b^ > a 2 > 0 , at r = 0 the first and third terms 
on the left hand side of (7.b) will be bounded, and only the 
second term will contribute to the stress singularity. In the 
second term a^ = 0, 0 < s < b^ , r > 0, and k 2 ^(r,s) is given 

by (8.d). Hence, the asymptotic value of the integral for small 
r may be obtained from ( 8 . d ) and (18. a, b). Using the definition 
(30) and ignoring the bounded terms, we then obtain 


MO) 


lim 

r+0 


/2 


“l 2y, 


T+E 


1 

IT 


2 

I 

k=l 


d 2k r 

(s+r) 1 


k-1 


■f 1 ( s ) ds 


= y 


* 

1 


9! (0) 


1 im /2 y,*r V, (r) , 
r+0 1 1 


l l " T+i^ 


2y i d 21 + ad 22 


si nwa 


m , P , 


( 3-2a^ ) ( 1 +m-j ) - (l“2ci|)(m-j+iC2) 
( 1 +m-| k-| ) (m^ +< 2 ) s i nira^ 


(31 .a, b) 


Similarly for b-j > a-j > 0, ^2 > a Z = we f inc * 


a. 


i o * g o ( o ) 

k ? ( 0 ) = lim /2 r x 1fi0 (r,Tr) = /I y ? 

2 r+0 160 2 /E7 


(32) 


where y 2 is obtained from (31. b) by interchanging the indices 
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1 and 2 in y j , m ^ , k . , (j = 1,2), and replacing by which, 
in turn, is obtained from (26. a) by substituting C£ ^ for c-^, 

(k = 1,2,3). It should be noted that the bielastic constant 

* * . 
y^ is identical to y given in [8] which was obtained by a 

different method. Needless to say, in this case for the remain 

ing three end points the expressions given by (28) and (29) 

will remain valid. 


For the third case in which the crack goes through the 

interface, i.e., for a-j = = 0, b^ >0, > 0, we define 

the normal and the shear components of the stress intensity 

★ 

factor as follows : 

k = lim r a x . Afl (r,Tr/2) , k = lim r a x. (r,ir/2) , 

e r ^. 0 jee r r + 0 jre 

(j = 1 or 2). (33) 

Thus, to obtain k Q and k r the expressions for the interface 
stresses are needed. These expressions may be obtained as 
definite integrals with the density functions f., (j = 1,2), 

J 


In this case, the asymptotic expressions of the stresses for 
small values of r and for 0 < 0 < it are of the form 

T. .(r,6) = h. .(e) + 0(r“ ), ( i , j = r,0; Re(a) > Re(-a 1 )) 

* J * J 

where h rtQ and h n are bounded and continuous in [0,tt], includ- 

t) t) r o 

ing at 0 = tt/2, and the bounded function h rr has a discontinu- 
ity at 0 = tt/2. Here, because of the physical importance of 
the interface and simplicity of the calculations, the stress 
intensity factors are defined in terms of the normal and 
shear stresses along the interface, i.e., k Q = k(O)h Q0 ( xr/ 2 ) , 

k r = k(0)h rQ (7r/2) . As the results of this section show, with- 
in a multiplicative constant, k ( 0 ) is nothinq but the constant 
g 1 ( 0 ) (which is related to g 2 ( 0 ) through (22)). 


22 



and the appropriate kernels which are evaluated from the origi- 
nal Mellin Transform solution (see [8]). After some manipula- 
tions, the components of the stress vector at the interface are 
found to be 


T iee^ r ’^ T 2ee( r ’J) i f £ "77 2 + 


b l A, s A 9 s(3r 2 -s 2 ) 


s +r 


( s 2 +r 2 ) 


2 — ]f 1 ( s)ds 


1 c. n Q ^ 


A~s A.s(3r 2 -s 2 ) 


s +r 


(s 2 +r 2 ) 


- — ]f 2 (s)ds 


_ n B,r B,r(r 2 -3s 2 ) 

T lre (r ’7 ) = T 2r0 (r »^ = ? / [ ~J~7 + ; — TT~ ] f ]( s ) ds 

0 s +r ( s 2 +r 2 ) 


1 f 2 r V , V< r2 - 3s2 > 

if ' L ~2 7 + 

0 s +r ( s 2 +r 2 ) 


2 — ]f 2 (s)ds » (34. a, b) 


A = - 

M 1 5 


A 2 " " l+m 1 K 1 ’ A 3 " 


y l 


m 2 +<^ 


A 4 " l+m 2 K 2 ’ 


^ +m i K 2 ^ y 2 ^ m i +lc 2^ 

^1 ( 1 +mi k-j ) (mi +<2 ) * ^2 ( 1 +m-j k-j ) (m-| +< 2 ) * 


y 1 (l+m 2 K 1 ) _ y 1 (m 2 +K 1 ) 

^3 ( 1 +m 2 <2 ) (ni2 + Ki ) * ^4 ( 1 ) (n^+K-i ) 


(35) 


The kernels in (34) can easily be expressed as the sum of terms 
containing 


1 1 
s+i r * s-i r * 


i r 

777? 


i r 


(s+ir) (s-ir) 


7 * 


(36) 


It is clear that since r > 0, at z = ± ir the functions G.(z) 

J 
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defined by (13) will be holomorphic and, following the same 
procedure which led to (18) in the previous section, the inte- 
gral s~in (34) may be evaluated in terms -of— G(z)- and -its deriv- 
ative (at z = ±ir). For example, 


, b l sf,(s)ds b l ' , 

¥ [ ~ 2 ~ 2 = 2? / ^ s+Tr + FTr ) f ]( s ) ds 


s +r 


• [ (e * 1cl/2 * e” 1 ^ 2 ) - 9 . V ( - 0) - L- + 


^1 r' 

b^ sirnra 


(r) 


9-] (0) 


1 


Bi a 

OK 1 cin™ r 

Zb^ sin^ - 


-4- + Mr) . 


/ ir ( 2 

o ( s+i r) 2 


1 

( s - i r ) 2 


)f ] (s)ds 



ag 1 (0) 


1 


^1 d ta r a 
1 Sln 2“ 


M + Mr) 


» 


(37. a, b) 


where and F 2 are bounded functions. Carrying out the neces- 
sary manipulations, from (33-37) we obtain 

■j 9 1 ( 0 ) 9 4 ( 8 ) 

ke = 1 "/if CAl ' (1 ’ 20l)A 2 :1 CA 3 - 

1 9 *1 ( 0 ) 9 o ( ® ) 

{ [b ’ - (, - 2 ° )B2] • if [B2 ■ (, - 2 ° )B < ]i> 

(38. a, b) 

where the constants g^(0) and g ^ ( 0 ) are linearly related through 
(22) and the constants A-j and B-j (i = 1,..,4) are given by (35). 


24 



5. 


CRACK SURFACE DISPLACEMENTS 


The crack surface displacements u-| 0 and u 20 are related 
to the density functions f-j and f 2 through (5). Hence, noting 
that Uj Q = 0 for r > bj , (j = 1,2), after evaluating f-| and f 2 
the displacements may be obtained from 


1 b] 

u l0 (r,TT-O) = j / f ] ( s )ds , 


l 2 

u 20 (r,+O) = - j j f 2 (s)ds 


(39. a, b) 


6. SOLUTION OF THE INTEGRAL EQUATIONS AND NUMERICAL RESULTS 

In deriving the integral equations (7), instead of the 
conditions (2.c) and (4.c) the conditions (6.a,b) are used; 
that is, the problem is formulated on the assumption that out- 
side the cuts (a.,b.), the displacement derivatives f i rather 
than the displacements Uj 0 , (j = 1,2), are zero. Thus, in the 
case of nonintersecting cracks (i.e., for a^ >0, a 2 > 0, or 
a^ = 0, a 2 > 0), the integral equations (7) must be solved 
under the si ngl e-val uedness conditions (6.c,d). Referring to 
[10], since the general solution of the system of singular 
integral equations [7] will contain two arbitrary constants,-^ 
the two additional conditions (6.c,d) are needed for a unique 
solution. On the other hand for the intersecting cracks, that 
is, for a^ = 0 = a 2 , there is only one si ngl e-val uedness con- 
dition, namely, -u i 0 ( 0 ,ir-0) = u 2q (0,+0), or 
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(40) 


U 1 u 2 

/ f, (s)ds = / f 2 (s)ds . 
o o 


The second condition which is necessary for a unique solution 
is provided by (22) relating g-|( 0 ) and g 2 ( 0 ) » or the functions 
f 1 and f 2 at r = 0 . 

For a^ = 82 = 0, even though the special cases of the 
resulting system of singular integral equations with generalized 
Cauchy kernels (7) have been considered in [13] and [14], there 
is no known method which can be used to regularize the integral 
equations of this type. To solve these equations we first nor- 
malize the intervals (a.,b.), (j = 1 , 2 ), through the following 

J J 

change in variables: 


2 s - (b,+a, ) 


2 s - ( b ? +a ? ) 
^2 = b 2 - a g 


( a j < s < b j ’ < < 1 ’ J = 1 .2) . 


(41 .a, b) 


Next we let 


fj(s) 


9j< s) 




a : 


(s-a)' j (b r s) 6j (H5 j ) aj (l-5 j ) 6j 


(42) 


giving 


9j(s) 


- ( 



a . + $ . 
J J 

) 




(43) 


Also defining 
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2r - (b.+a.) 

t = sL - J. 


b . - a . 
J J 


W - P j < r ) , 


( a j < r < b j * - 1 < tj < 1 » j = 1 .2) , 


( 44 ) 


and suppressing the subscripts in £ . and t., the system of 

J J 

integral equations (7) may be expressed as 


1 

IT 


| C + ^ P,(t) , 


(-1 < t < 1), (45) 

where 

Wj(£) = ( 1 +£) aj (l-£) , (j = 1,2) (46) 

are the fundamental functions of the system, and the kernels 
K.j j may be obtained from k ^ j by means of appropriate transfor- 
mations. Now, observing that the fundamental functions w.(£) 

J 

of the system of integral equations are weights of the Jacobi 

( - B . , - a . ) 

polynomials P n J J (£), the equations may be solved numer- 
ically for the unknown functions <}>• by using the Gauss-Jacobi 

J 

integration formulas described in [15]. After evaluating 4> • , 

J 

the values of g . at the end points, which are needed to calcu- 
J 

late the stress intensity factors, may be obtained from (43) 

i n terms of <J> . (±1 ) . 

J 

For the material pairs Aluminum-Epoxy some of the calcu- 
lated results for the plane strain case are shown in Tables 
3-5 and Figures 4-13. For b-j 0 or b 2 -»■ 0 (i.e., for the case 
of a single crack terminating at the interface), the results 
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of [8] are obtained as the limiting values. In the numerical 
examples given in this paper it is assumed that the composite 
m e d i u m — is - lo a ded - pa r a VVe~l~t o — t h e — i n t e r face” and” awayf rom^the 
crack region in such a way that the stresses p-|(r) and PgCr) 
are approximately constant. Thus, in the singular part of the 
problem considered in this paper, the crack surface tractions 
will be constant pressures -p-| and -p 2 which are related 
through ( 1 ) . 

In the example shown in Table 3 and Figures 4 and 5, the 
material 1 is aluminum, material 2 is epoxy, one of the crack 
tips b-| is fixed at a unit distance from the interface, and 
b 2 varies between 0 and 2 b ^ . In Figure 4 as well as 6 and 8 
the constants k.| and k 2 refer to k(b-j) and k^), respectively. 

Table 3. The Stress Intensity Factors for Through Crack, 
a^ = a 2 = 0. Material 1: Aluminum, Material 2: 

Epoxy (vj = 0.3, V 2 = 0.35, E i / E 2 = 22.447); 
b 1 = 1 = constant, i = (b-|+b 2 )/ 2 , b 2 variable, 

.(8} = 8 2 = °- 5 » « = 0.273692). 


b 2 /b, 

k(b ] )/ P] /T 

k(b 2 )/p 2 /Z 

k e /p,n“ 

k r /p } Z a 

ro 

1 .3552* 

00 

-h - OO 

00 


1.40374 

4.36065 

-0.22477 

0.08215 

ns 

1.33240 

2.13917 

-0.08113 

0.03184 

H9 

1.23789 

1 .5561 5 

-0.05906 

0.02158 

0.75 

1 .1 5887 

1.31166 

-0.04636 

0.01694 

1.00 

1.09312 - 

1.17874 

-0.03839 

0.01403 

1.25 

1.03767 

1 .09693 

-0.03266 

0.01194 

1 .50 

0.99020 

1 .04261 

-0.02822 

0.01032 

1.75 

0.94902 

1.00462 

-0.02462 

0.00900 

2.00 

0.91286 

0.97701 

-0.02160 

0.00789 


*Taken from Reference [8], 
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Table 4. 


The Stress Intensity Factors for Through Crack, 
a-j = &2 = 0- Material 1: Aluminum, Material 2: 
Epoxy, 2Z = b i +b 2 = 2 = constant, c = (b 2 -b-j )/2 
variable, ( 6 -j = = 0.5, a = 0.273692). 


c/l 


' k(b 2 )/P 2 /^ 

k e /p i a0t 

Vp ^ 01 

■m 

1 .3552* 

-y 00 

-> -00 

-y 00 


1 .40333 

5.69903 

-0.34124 

0.12472 


1.37532 

2.76781 

-0.11960 

0.04371 

-0.50 

1 .29929 

1 .87137 

-0.07423 

0.02713 

-0.25 

1.20448 

1 .43765 

-0.05310 

0.01941 

0.00 

1 .09312 

1.17874 

-0.03839 

0.01403 

0.25 

0.96213 

1 .01590 

-0.02575 

0.00941 

0.50 

0.80281 

0.91834 

-0.01293 

0.00473 

0.75 

0.59068 

0.87340 

0.00313 

-0.00114 

0.95 

0.29854 

0.87519 

0.02903 

-0.01061 

1.00 

-*■ 0 

0.8827* 

00 

-> -00 

*Taken 1 

From Reference [ 8 ]. 


Table 5. The Stress Intensity Factors for Through Crack, 
a-j = a£ = 0. Material 1: Epoxy, Material 2: 
Aluminum; b-| = 1 = constant, Jl = ( b-j +bg ) / 2 , 
b 2 variable, (a = 0.273692 , 6 -j = 62 = 0.5). 


b 2 /b l 

k ( b 1 )/p 1 /l 


k 0 / p l«. a 


m 

0.8827* 

-► 0 

-y 00 

-y 00 

■89 

0.87141 

0.39001 

0.44554 

0.16283 

mm 

0.89435 

0.72689 

-0.16129 

-0.05895 

■89 

0.97701 

0.91286 

-0.48485 

-0.17720 

0.75 

1.07647 

1 .02106 

-0.69736 

-0.25487 

1 .00 

1.17874 

1.09312 

-0.86170 

-0.31493 

1.25 

1 .27901 

1.14476 

-0.99895 

-0.36509 

1.50 

1 .37559 

1 .18357 

-1.11885 

1 

-0.40891 

1.75 

i 1 .46797 

1 .21377 

1 -1.22669 

-0.44832 

2.00 

1.55615 

1.23789 

-1.32567 

-0.48449 

*Taken 1 

From Reference [8]. 
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Figure 4. Stress intensity factors for a crack going through 
the interface (materials: 1 aluminum, 2 epoxy; 
b, = 1 = constant, * = (bj+bg)^ variable). 



Table 4 and Figures 6 and 7 show the results for the same mate- 
rial combination. In this second example it is assumed that 
the total length 21 = b^+ b 2 of the crack is constant, and the 
distance c = (bg-b^ )/2 from its center to the interface is 
variable. The stress intensity factors for the case of epoxy- 
aluminum are shown in Table 5 and Figures 8 and 9. 

The tables also show the limiting values of the stress 
intensity factors for b 2 -*■ 0 and c -*■ ±1 . These limiting values 
are obtained by observing that in the case of a through crack 
(i.e., for b. > a. =0, (j = 1,2)), for the material pair 

J J 

aluminum-epoxy (21) gives the power of the stress singularity 
as a = 0.2737. On the other hand, for the limiting case of a 
crack tip terminating at the interface, i.e., for b^ = 0 or 
bg = 0, and for the material pair under consideration, from 
(26. a) the power of the stress singularity is calculated to 
be: 

= 0.8258 , crack in aluminum, 

ot.| = 0.3381 , crack in epoxy. (47. a, b) 

Thus, for example, in the example shown in Table 3 and Figures 
4 and 5, as b 2 -*• 0 the powers of singularity become a -*■ a-j , 

$2 ■* a-j • Therefore, for this limiting case, since = 0.8258 
> a = 0.2737 and a-j > = 0-5, at r = 0 the stresses would 

have a stronger singularity and consequently the stress inten- 
sity factors k Q and k . which are defined on the basis of r" a , 
and k ( b g ) , which is defined on the basis of r"®‘®, would all 
go to infinity. Similar arguments apply to the limiting values 
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Figure 6. Stress intensity factors for through crack 
(materials: 1 aluminum, 2 epoxy; 21 - 2 - 
constant, c variable). 
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Figure 7. Stress intensity factors for normal and shear 
stresses on the interface for through crack 
(materials: 1 aluminum, 2 epoxy; 2l = 2 = 
constant, c variable). 



0.50 1.00 1.50 

Jt 

Figure 8. Stress intensity factors for through crack 
(materials: 1 epoxy, 2 aluminum; b-j = 1 = 
constant, l variable). 
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Figure 9. 


□ 


X 

Stress intensity factors for normal and shear 
stresses on the interface for through crack 
(materials: 1 epoxy, 2 aluminum; b-j = 1 = 
constant, X variable). 


shown in Tables 4 and 5. 

From Tables 3-5 and Figures 5, 7, 9 it is seen that the 
normal component of the stress intensity factor, k 0 , (hence, 
the normal component of the bonding stress along the interface 
in the neighborhood of r = 0) may be positive or negative de- 
pending on the relative penetration depths b-j and b 2 of the 
crack into the adjacent materials. As the ratio bg/b^ in- 
creases, k 0 increases if p-j > v> 2 (Figures 5, 7), and decreases 
if P 1 < y 2 ( Figure 9) . 

Figure 10 shows a sample distribution of the density 
functions f-| and f 2 defined by 

f-j(s)-F-|(5-j), (0 < s < b-j , -1 < Ci < 1), 

f 2 ( s ) = f 2 (C 2 ) » (0 < s < b 2 , -1 < C 2 < 1)- (48. a, b) 

The calculated values shown in the figure correspond to the 
example given in Table 3 with bg/b^ = 0.05 (i.e., material 1 
is aluminum, material 2 is epoxy). The important feature of 
the distribution of the density functions fi and f 2 is that at 
r = 0, because of the common r" a type singularity, they become 
infinite, -fi and f 2 having the same sign (i.e., at r = 0 the 
derivative of the crack surface displacement is discontinuous 
and unbounded) . 

Sample results for the crack surface displacement calcu- 
lated from (|39) are shown in Figures 11-13, where the variables 

! 

are again the normalized quantities C] and C 2 defined by (41). 
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Figure 10. Distribution of the density functions fi(r) = 
Fi ( Ci ) » "f o ( r ) = f 2^2^ (materials: 1 aluminum 
2 epoxy; Dg/b-j = 0.05; Cj = ( 2r-b j -a j )/ ( bj -a j ) 


The results shown in Figures 11, 12, and 13 correspond to the 
examples given in the Tables 3, 4, and 5, respectively. In 
Figures 11 and 13 the length unit for displacements is b-j and 
in Figure 12 it is l = (b-j+bg)/^. A careful examination of the 
crack surface displacements in the adjacent materials around 
r * 0 would indicate that the crack opening in the medium with 
the higher modulus is greater than that with the lower modulus, 
both having an infinite slope. 

The results given in this paper may be used in the analy- 
sis of fracture propagation in the composite medium by using 
any of the standard fracture criteria if the crack extension 

takes place at the singular point r = b^ or r = where the 

- 1/2 

singularity is of the form (r-b.) ' . However, in studying 

J 

a possible fracture initiation at r = 0 (where the singularity 
is of the form r” a ) , i.e., initiation of a debonding crack 
along the interface or a cleavage crack in one of the adjacent 
materials, a criterion similar to that outlined in [8] may 
have to be used. 
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Figure 13. Crack surface displacements for through crack 

(materials: 1 epoxy, 2 aluminum; b, ■ 1 = constant) 
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